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摘要 :以 鞍 形 膜 结构 为 研究 对 象 ,研究 了 冲击 载荷 作用 下 鞍 形 膜 结构 的 振动 问题 。 首 先 推导 了 冲 
击 载荷 作用 下 的 结构 振动 非 线 性 微分 方程 ,该 方程 中 考虑 了 膜 材 的 正 交 异性 ,空气 阻尼 以 及 膜 材 变 
形 的 几何 非 线 性 。 然 后 假设 一 小 球 以 茶 确 定 的 冲击 速度 差 冲击 膜 面 ,模拟 冲击 荷载 。 最 后 利用 随 
机 摄 动 法 求解 方程 ,得 到 鞍 形 膜 结 构 振 动 响应 的 位 移 函 数 ,通过 该 函数 分 别 对 膜 面 预 张力 、 拱 跨 比 
上 以 及 小 球速 度 差 进 行 参 数 分 析 。 同 时 ,对 鞍 形 膜 结 构 的 振动 进行 数值 模拟 ,并 与 解析 理论 作对 比分 
< 析 。 结 果 表 明 鞍 形 膜 结 构 振 动 位 移 会 随 着 膜 面 预 张力 或 拱 跨 比 的 增 大 而 减 小 ,同时 也 会 随 着 速度 
) 差 的 增 大 而 增 大 。 由 此 可 知 该 理论 模型 研究 鞍 形 膜 结构 在 冲击 载荷 作用 下 的 振动 特性 是 可 行 的 ， 
得 到 的 结果 能 够 准确 地 预测 膜 结 构 响 应 规律 。 
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Abstract : Taking the saddle membrane structure as the research object,the vibration of the saddle mem- 
brane structure under the impact load is studied. First, the nonlinear differential equation of structural vi- 
bration under impact load is derived ,in which the orthotropic of membrane material , air damping , and geo- 
metric nonlinearity of membrane deformation are considered. Then,it is assumed that a small ball impacts 
the membrane surface with a certain impact velocity difference to simulate the impact load. Finally , the sto- 
chastic perturbation method is used to solve the equation , and the displacement function of the vibration re- 
sponse of the saddle membrane structure is obtained. At the same time ,the vibration of the saddle mem- 
brane structure is numerically simulated and compared with the analytical theory. The results show that the 


vibration displacement of saddle membrane structure decreases with the increase of membrane pretension 
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or arch span ratio ,and increases with the increase of velocity difference. It can be seen that the theoretical 


model is feasible in studying the vibration characteristics of saddle membrane structure under impact load, 


and the results can accurately predict the response law of membrane structure. 
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膜 材 由 于 其 质量 轻 、 柔 韧性 好 , 现 被 广泛 应 用 于 
大 跨度 空间 建筑 结构 当中 ,并 且 工程 师 也 常用 正 交 
各 向 异性 织物 膜 材 设计 成 空间 曲面 形状 ( 伞 形 和 较 
形 ) 的 屋顶 或 窒 顶 结构 。 因 为 刚度 低 、 质 量 小 的 特 
点 , 膜 结构 在 受到 外 界 风雨 冰雹 和 风 卷 产物 等 荷载 
作用 时 ,会 十 分 敏感 , 易 产生 较 大 的 振动 变形 ,甚至 
结构 失效 。 而 实际 环境 中 膜 结构 受 风 荷载 作用 更 为 
广泛 , 且 当 今 比 较 典 型 膜 结构 工程 事故 中 ,有 多 起 是 
爱 风 荷载 作用 引起 的 ,因此 许多 学 者 对 风 荷 载 问题 
进行 了 较为 广泛 的 研究 。 而 膜 结构 的 冲击 问题 就 相 
对 蝎 得 不 那么 突出 ,相关 的 研究 也 较 少 ,但 是 膜 结构 
了 竹 的 设计 、 施 工 与 维护 , 仅 考虑 风 荷 载 作 用 是 不 全 
三 = 处 合理 的 ,如 较 大 的 冰 夫 冲击 会 使 膜 结构 产生 较 
龙岗 振动 变形 以 及 应 力 松弛 "1。 因 此 ,为 了 充分 考 
VECES PIE BOUE AL ,保证 膜 结构 工程 的 可 靠 性 , 研 
究 膜 材 在 冲击 载荷 作用 下 的 振动 问题 是 十 分 必 
= 目前 , 膜 结构 相关 问题 的 研究 较为 广泛 ,包括 各 
种 形状 模型 的 研究 ,以 及 确定 性 的 与 非 确定 性 的 力 
学 研究 等 。 张 华 等 ”使 用 动力 松弛 法 对 薄膜 结构 进 
行 静 载 分 析 , 提 出 一 种 方法 处 理 争 折 单 元 以 确保 荷 
载 分 析 的 可 靠 性 。LIU 等 ”考虑 了 薄膜 几何 非 线 
性 6 研究 了 正 交 各 向 异性 膜 薄 膜 阻尼 作用 以 及 冲击 
载荷 作用 等 问题 ,得 到 了 薄膜 振动 频率 . 振 型 以 及 位 
移 解析 ,为 平面 膜 结构 设计 、 膜 面 预 张力 测量 问题 提 
供 了 理论 依据 。ZHENG 等 ”通过 试验 得 到 了 冲击 
载荷 试验 的 统计 特性 ,通过 随机 振动 理论 研究 了 膜 
材 的 随机 振动 问题 ,得 到 了 振动 位 移 响应 均值 函数 
的 计算 方法 ,为 随机 振动 问题 的 进一步 研究 奠定 了 
基础 。LI 等 “基于 随机 脉冲 理论 ,获得 了 随机 脉 
冲 风 荷 载 特性 ,研究 了 正 交 各 向 异性 平面 膜 在 随机 
脉冲 风 作 用 下 的 随机 振动 问题 ,为 随机 振动 控制 和 
随机 动 载荷 作用 下 膜 结 构 的 设计 提供 了 理论 框架 。 
郭 学 敏 等 思 针 对 工程 中 电子 器 件 的 膜 结构 热 应 力 与 
热 应 变 会 严重 影响 器 件 结构 与 功能 的 问题 进行 了 研 
究 , 解 决 了 困扰 有 限 元 分 析 的 奇异 点 问题 ,结果 对 恶 
劣 工作 环境 下 电子 器 件 的 薄膜 设计 具有 一 定 的 借鉴 


意义 ,为 膜 结 构 热 力学 研究 提供 了 参考 。 对 空间 曲 
面 模型 的 膜 结构 ,由 于 微分 控制 方程 中 曲率 的 影响 ， 
大 大 增加 了 数学 上 的 复杂 性 ,很 难 求 得 曲面 膜 振 动 
问题 的 精确 解析 解 。 因 此 ,多 采用 近似 方法 求解 近 
似 的 解析 解 , 或 采用 试验 和 数值 模拟 方法 进 研 究 。 
如 XU 等 中 LIU 等 "通过 建立 基于 冯 “' 卡门 大 挠 
度 理论 和 达 朗 贝尔 原理 的 风 振 耦合 控制 方程 ,研究 
了 正 交 各 向 异性 鞍 形 膜 结构 的 非 线 性 风 振 气 动 稳定 
性 。 通 过 判断 系统 特征 方程 的 稳定 性 ,确定 了 发 散 
失 稳 临 界 风速 。 孙 芳 锦 等 … 研究 了 风 与 柔性 结构 
流 固 耦 合作 用 的 强 耦 合 整 体 方程 ,采用 预 处 理 方法 
对 经 典 二 维 与 三 维 流 固 耦 合 问题 进行 了 计算 分 析 ， 
结果 表明 该 方法 可 使 强 耦 合 整体 方程 的 求解 在 计算 
精度 和 效率 上 得 到 较 大 提升 ,证 明 该 方法 适用 于 风 
与 柔性 结构 的 流 固 耦 合 分 析 ,为 风 与 膜 结 构 的 耦合 
作用 研究 提供 了 方法 。 

综 上 ,作为 膜 结构 常用 形式 之 一 的 鞍 形 膜 结构 ， 
其 冲击 问题 还 鲜 见 相关 报道 。 因 此 ,本 研究 以 鞍 形 
膜 结构 为 研究 对 象 ,利用 求解 各 种 非 线 性 微分 方程 
的 强 有 力 方法 一 一 随机 摄 动 法 ,求解 得 到 其 在 冲击 
载荷 作用 下 的 振动 理论 解析 ,通过 该 解析 得 到 鞍 形 
膜 结构 的 响应 位 移 函 数 ,并 分 析 膜 结构 相关 参数 对 
鞍 形 膜 结构 的 动力 响应 的 影响 规律 ,通过 Ansys/Ls- 
Dyna 有 限 元 软件 进行 数值 模拟 ,与 理论 作对 比分 
析 。 最 终 对 比 结果 表明 理论 解析 和 数值 结果 吻合 较 
好 ,能 够 准确 地 预测 鞍 形 膜 结构 在 冲击 载荷 作用 下 
的 振动 特性 ,得 到 冲击 荷载 作用 下 的 膜 结构 响应 规 
律 ,通过 该 规律 可 为 确定 膜 结构 相关 参数 的 设计 ( 如 
拱 跨 比 选择 、 预 张力 确定 与 测量 ) 和 膜 结构 施工 与 维 
护 提供 参考 ,也 为 鞍 形 膜 结构 的 随机 振动 与 可 靠 度 
分 析 提 供 理论 基础 。 


1 边界 条 件 与 控制 方程 


假设 鞍 形 膜 结 构 的 四 条 边 为 简 支 支承 。 以 膜 表 
面 的 一 个 角 点 o 为 原点 建立 空间 直角 坐标 系 ,x\y 方 
向 尺寸 分 别 为 a、b;x、y 方向 预 张力 分 别 为 No, M 
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No, 5A 点 为 膜 面 中 心 点 ,B、C、D 点 分 别 膜 面 各 位 置 
测 点 汰 状 分 别 是 结构 xy 方向 的 跨 中 拱 ;假设 一 小 
球 卫 以 初速 度 vy 冲击 膜 面 中 心 点 4, 后 以 速度 v BE 
反射 出 膜 面 , 如 图 1 所 示 。 


Ce 。 图 1 四 边 简 支 正 交 异 性 鞍 形 腊 结构 模型 
LC) Fig.1 Quadrilateral simply supported orthotropic 


e : saddle membrane structure model 
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应 力 边 界 条 件 为 
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其 中 : M t) 为 横向 冲击 产生 的 位 移 响应 ; o (o, 
7 站 为 薄膜 应 力 函 数 六 为 任意 时 间 。 

较 形 膜 结 构 在 冲击 荷载 作用 下 的 偏 微 分 控制 方 
程 根据 参考 文献 ”" ,得 到 
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式 中 :h 为 腊 材 厚度 ;C 为 阻尼 系数 ;p 为 膜 材 面 密 
BE SERIE, 4319178 x I y 向 弹性 模 量 ;Ko, 和 Ko, 分别 
为 x 和 y 向 的 曲率 。 

假设 在 一 定时 间 上 和 一 定位 置 (xo ,yo ) 对 膜 面 施 
加 一 个 振幅 为 4 的 冲击 载荷 ,表示 为 

P(x,y,) = AB(0)8(x -,y - y) (6) 

式 中 :6 为 犹 拉克 脉冲 函数 ;4 为 冲击 荷载 幅 值 。 


2 求解 非 线 性 微分 方程 


根据 板 沈 振动 理论 
式 (2) 的 位 移 函 数 为 


w(x,y,t) = Y Y sn (79) (7). 了， (t) 


mzlnzl 


,假设 满足 位 移 边 界 条 件 


= DD Wp) Tl) (7) 


n=1 nzl 


RP: W, Gy) 为 振 型 函数 ; T, C0) 为 未 知 函数 ;m 
和 nn 为 整数 ,分 别 表示 x 和 y 方向 上 的 正弦 半 波 数 。 

将 式 (7) 代 入 相 容 方程 式 (5 ) ,再 结合 边界 条 件 
式 (3) ,可 以 得 到 应 力 函 数 为 
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TEX C7) ~ (8) 代 入 弹性 曲面 微分 方程 式 (4) , 
式 中 7(1t) 用 7 表示 , 则 运用 件 辽 金 法 解 看 可 得 
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上 式 引 入 小 参数 e = n^ 7ab 并 化 简 为 非 线 性 振 
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动 微分 方程 得 (9) 
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205 随机 摄 动 法 求解 
AIRO) 是 关于 TO) 的 非 线 性 振动 微分 方程 ,本 
三 完 利用 随机 振动 理论 中 的 随机 摄 动 法 进行 求解 。 
由 宁 膜 的 厚度 远 远 比 边 长 小 , 即 e = h/ab <1 ,所 
VIX TO) 展开 成 关于 小 参数 e 的 宕 级 数 ,得 到 
S TORT +eT e eT, +oe) (10) 
COR C10) 代入 式 (9) ,并 令 e 的 同 次 宕 各 项 系 
BETE ,可 以 得 到 一 组 线性 微分 方程 , 即 
35-287, 4 uiT. o Xp 
T, -2BT, +o T, =- auo T. 
T, + 2BT, + oT, 2 - 2a o», T,T, - aso T, 


(11) 


HRL) nl, TCO) 的 每 一 个 展开 式 都 满足 一 
个 输入 的 线性 微分 方程 。 根 据 平稳 随机 过 程 理论 ， 
可 以 对 其 进行 杜 哈 梅 积 分 求解 , 即 


T= | hG)XG - ds, 


T, = | hegi G - sds, 


T, =Í hs g(t - s)ds (12) 


式 中 : h(s) 为 脉冲 响应 函数 ; e CO = - aioo70 s 


2 258 
g2 7 -2a,e ToT, - aw To o 


一 个 样本 ,因此 可 以 是 一 种 确定 性 非 线 性 振动 问题 。 
当 小 球 撞击 膜 表 面 时 , 膜 表 面 与 球 之 间 的 碰撞 
是 一 对 作用 力 和 反作用 力 。 假 设 球 的 入 射 速 度 为 
vo ,反射 速度 为 w 。 根 据 D'Alembert's 定理 ,可 以 确 
定 膜 表 面 的 冲击 力 为 
MDC = Mv, - Mv, >A = M8(v) (13) 


由 式 (9) 可 得 到 作用 于 贰 形 膜 的 冲击 载荷 过 程 


Xi) E Jus (2*7 )sin(* 770)» 
a 


(14) 
根据 文献 15], 由 式 (11) 推 导 脉 冲 响应 函数 
h(s) 为 


e 
h = —c———53in( 3g (15) 
PS cg Weser 

2.2 位 移 响 应 函数 


关于 T(t) 的 非 线性 振动 微分 方程 , 随 着 e 的 阶 
数 越 大 运算 越 复杂 ,实践 中 党 只 计算 到 s WRR, 
而 高 次 摄 动 也 只 能 改善 小 e 的 近似 程度 ,对 于 大 e 
时 高 次 摄 动 可 能 使 近似 程度 变 坏 : 鸣 。 因 此 将 式 
(14) ~ (15) 代 入 式 (12) 的 第 一 个 积分 项 可 得 
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T, = 。 sin[ Z^ Jsin[ T. 
pab Vos - B. 


sin(V/ — gii) (16) 
再 将 式 (16) (12) 中 的 gi G) 代入 式 (12) 的 第 
二 个 积分 项 可 得 


- 16a,0, M^8 G)* sin? (T9 sin (270) 
a 


b 
Jäs -g Coab?? (wo -p» 


T, = 


B Vo - B'sin( Joz gu) e - em 
36^ - ABoj i 
Va - p'e” PauB -oge 2* 
2cos( Jas -BU 
TBE -46m 
(17) 


将 式 (10) ,(16) (17) 代 入 式 (7) , RRE 
腊 结 构 在 膜 表面 任意 点 处 位 移 的 近似 解 为 


© a g 
wey, D= 之 È sin( 77 sin( 2P) T, +e- T) 
Se. m-lnzl a 


(18) 
AI 将 鞍 形 膜 结构 的 初始 表面 函数 式 (1) 与 式 (18) 
进 符 半 加 ,可 以 得 到 冲击 荷载 作用 下 进 形 膜 结构 的 
曲面 函数 。 


y C QD! fi Gr - (O/2)* 
dd (a/2? qv 


(19) 
-三 式 (18) ~ (19) rh à Co) 是 小 球 冲击 膜 面 过 程 的 
速 语 差 。 取 确定 性 的 速度 差 8(v) 输 入 ,可 以 得 到 确 
定性 的 动态 响应 输出 。 


3 ABAH 


以 工程 中 常用 膜 材 为 例 ,具体 参数 取 :a =b = 
I m,E, =1 400 MPa, E, =900 MPa,p =1.7 kg/m’, 
h 20.001 m, HJE C 2120 N * sm, F EXA 
JE BTE S E TE RS) np Rom EHI P IREE o 


3.1 位 移 时 程 曲线 

BUBt E f, - f, =0.1, 预 张力 No = No = 
10 kN/m, H 3X (2) 可 得 , Ky = 0. 8 m^, Ko = 
-0.8 m ,假设 质量 为 M=25 x10 7 kg 的 小 球 冲 
击 膜 面 中 心 点 4。 根据 式 (18) ,冲击 点 4 的 前 三 阶 
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位 移 以 及 前 三 阶 位 移 钱 加 时 程 曲 线 如 图 2 所 示 。 


000 00 004 006 008 010 
时 间 /s 


图 2 冲击 点 4 前 三 阶 位 移 时 程 曲 线 
Fig.2 The first three-order displacement 


time-history curve of impact point A 

从 图 2 可 以 得 出 以 下 结论 。 

1) 最 大 位 移 随 着 振动 阶 数 的 增加 而 逐渐 减 小 ; 
而 振动 周期 也 随 着 振动 阶 数 的 增加 而 减 小 。 

2) 随 着 时 间 的 增加 ,由 于 阻尼 的 存在 , 单 阶 振 
幅 逐 渐 减 小 ,直至 为 零 。 

3) 总 位 移 振动 波形 更 接近 实际 的 不 规则 性 , 同 
时 呈现 出 各 阶 振动 波形 周期 振幅、 衰减 等 特点 。 且 
当 小 球 冲击 膜 面 后 ,4 点 做 有 阻尼 自由 振动 。 

同时 对 膜 面 另 外 3 个 测 点 B.C、D 分别 进行 了 
前 3 阶 模 态 三 加 ,得 到 总 位 移 的 时 程 曲 线 ,并 与 4 点 
作 比 较 , 如 图 3 所 示 。 


位 移 /mm 


RH s : 2 
0.00 0.02 0.04 0.06 0.08 0.10 
时 间 /s 


图 3 各 测 点 前 三 阶 位 移 时 程 曲线 
Fig.3 The first three-order displacement time-history 
curve of each measuring point 

从 图 3 可 以 看 出 ,小 球 冲击 4 点 后 , 膜 表面 产生 
振动 变形 。 而 B.C、D 点 产生 的 振幅 均 比 4 点 小 。 
B.C 两 点 由 于 处 在 位 置 的 曲率 不 同 ,两 点 各 自 的 振 
幅 以 及 周期 有 一 定 差异 。D 点 是 由 正 向 的 小 振幅 开 
始 的 有 阻尼 自由 振动 ,这 与 其 他 三 点 有 所 差别 ,这 是 
由 于 振动 波形 传递 时 , 波 的 干涉 因 各 向 异性 以 及 预 
张力 的 不 同 导致 D 点 为 正 向 小 振幅 。 


(a) 一 阶 £-0.002 s 


19 99 


(c) 三 阶 10.002 s 
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2 v 
10 0.0 


(b) 二 阶 1=0.002 s 


S, 
10 9.0 


(d) Bi =Æ) (0.002 s 


图 4 f IMS Bi ZNE d 


Fig.4 Vibration modes of each order and the first three superposition modes 
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30. EE BE ERES 


.全 将 材料 和 几何 参数 代入 式 (19) ,为 获得 较为 明 
TIRE , 取 小 球 质量 为 M =25 x10 7 kg。 可 以 得 
到 交 阶 的 振 型 ,也 可 以 得 到 各 阶 的 又 加 振 型 。 振 型 
如 图 4 所 示 , 华 标尺 寸 为 m. 

加 从 以 上 振 型 结果 可 以 看 出 ,利用 鞍 形 膜 曲 面 函 
数 式 (19) 可 以 方便 地 计算 出 寺 形 膜 结构 在 冲击 载 
荷 作用 下 的 各 阶 振 型 以 及 总 释 加 振 型 。 振 动 幅 值 随 
振动 阶 数 的 增加 而 减 小 , 即 低 阶 振 型 的 影响 大 于 高 
阶 振 型 。 


4 数值 模拟 


壕 形 膜 结构 在 冲击 载荷 作用 下 的 动力 响应 数值 
分 析 是 典型 的 预 张力 结构 瞬 态 分 析 问 题 。 针 对 这 类 
问题 ,本 研究 利用 了 通用 显示 动力 学 有 限 元 分 析 软 
件 ANSYS/LS-DYNA 对 冲击 载荷 作用 下 贰 形 膜 结构 
振动 进行 了 模拟 " ,并 与 理论 结果 进行 对 比分 析 。 

求解 方法 采用 隐 式 - 显 式 连 续 解 。 在 隐 式 计算 
中 , 膜 采 用 了 Shell181 单元 ,小 球 采 用 了 Solid185 单 


元 。 而 在 显 式 动态 分 析 时 ,将 隐 式 分 析 中 的 相应 单 
元 转换 为 Shell163 单元 、Solid164 单元 。 网 格 划 分 
中 ,采用 了 映射 三 角形 单元 来 生成 膜 表 面 , 映 射 六 面 
体 实体 单元 来 生成 小 球 的 网 格 , 冲 击 点 膜 面 网 格 大 
小 为 0.25 em ,小 球 网 格 大 小 为 1 em, 网 格 划 分 效果 
如 图 5。 而 数值 模拟 与 理论 取 相 同 材料 参数 ， 
见 表 1。 


图 5 小 球 单元 与 膜 单元 网 格 划分 
Fig.5 Mesh division of small ball element 


and membrane element 
对 于 小 球 冲 击 膜 面 的 速度 差 值 的 模拟 ,通过 在 
显示 分 析 关 键 字 选项 中 ,设置 Solid164 单元 Z 方向 
的 不 同 初速 度 , 进 行 LS-DYNA 求解 ,然后 在 LS-Pre- 
Post 后 处 理 程序 中 查看 Solid164 单元 的 速度 时 程 曲 
线 , 得 到 不 同 初速 度 下 的 速度 差 值 。 这 样 ,设置 不 同 
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的 初速 度 以 满足 速度 差 , 达 到 对 小 球 冲 击 膜 面 速度 


差 值 模拟 。 因 此 ,不 同 速度 差 所 对 应 小 球 的 模拟 初 
速度 以 及 弹出 膜 面 的 反弹 速度 见 表 2。 


表 1 材料 属性 
Tab.1 The material properties 
B 厚度 密度 p/ ”弹性 模 量 ”弹性 模 量 ” 泊 松 比 
单元 类 型 
h/m (kg-m 3)  E,/MPa E,/MPa NUXY 
SHELLIAI 0.001 1 700 1 400 900 0.34 
SOLIDI85 一 900 53 53 0.3 


R2 ， 同 速度 差 所 对 应 的 模拟 速度 
Tab.2 Simulation speed corresponding to 


different velocity differences 


速度 差 /(m * s7!) 


速度 
38 39 40.7 4 43 


XPBLEDEEHE/(m- s^!) 35.1 36.7 37.18 238.7 39.5 


HEEE (m-s!) 3.0 3.2 36 3.9 3.9 


SCR S HOO LEA 


CS BOE RES fe nti Bo MEA F K T fe RUE 
FERE f£ 71 ET SUBE rf ^ "E B ACRI. A 
此 > 位 移 被 视 为 振动 分 析 和 动态 设计 的 重要 参数 。 
大 暗 详细 讨论 了 预 张 力 、, 拱 跨 比 ,冲击 载荷 速度 差 等 
变量 对 位 移 结果 的 影响 。 参 数 取 a - b = 1 nm， 
E31 400 MPa, E, = 900 MPa, p = 1. 7 kg/m', 
M e. 001 m,C - 120 N - s/m,M 0.25 x 10? kg。 


SE 预 张力 对 位 移 最 大 值 的 影响 


号 选取 预 张力 分 别 为 5.6、7、8.9、10 kN/m; fis 
比 与 速度 差 均值 取 f. =f, 20.1, 8(v) = 40.7 m/s; 
计算 出 所 选 预 张 力 的 鞍 形 膜 在 振动 时 产生 的 最 大 位 
移 , 结 果 见 表 3 和 图 6。 

表 3 不 同 预 张力 下 的 最 大 位 移 比较 

Tab.3 Comparison of the value of maximum displacement 

under different pretension 
预 张力 /kN 
5 6 7 8 9 10 


d: H 
25H 


p 


理论 结果 /mm 9.94 9.64 9.35 9.07 8.79 8.52 
数值 结果 /mm 9.35 9.11 8.86 8.74 8.65 8.39 
相对 差 /% 5.94 5.50 5.24 3.42 1.59 1.56 


从 表 3 和 图 6 可 以 看 出 , 随 着 预 张力 的 增加 ,位 
移 值 呈 下 降 趋 势 , 预 张力 从 5 kN/m 到 10 kN/m 之 
间 , 膜 面 位 移 理论 解 减少 14% , 而 数值 解 减少 10% 。 


;1 这 表明 预 张力 的 变化 对 鞍 形 膜 结 构 抵 抗 振动 变形 有 


明显 的 影响 ,能 为 膜 面 提供 一 定 抵抗 变形 的 刚度 ,所 
以 这 对 膜 结构 设计 来 说 是 一 项 重要 参数 。 同 时 , 理 
论 结果 和 数值 结果 之 间 的 相对 差 逐渐 降低 ,最 大 相 
对 差 达 5.94% ,表明 预 张力 变化 对 鞍 形 膜 结构 振动 
的 影响 ,在 数值 解 与 理论 解 上 的 变化 是 吻合 的 。 


10.0 


一 一 理论 解 


一 * 一 数值 解 


^5 6 7 8 9 10 
预 张 力 /kN 
图 6 不 同 预 张力 下 的 最 大 位 移 理论 解 与 数值 解 对 比 


Fig.6 Comparison between theoretical solution and 


numerical solution of maximum displacement 


under different pretension 
5.2. 拱 跨 比 对 最 大 位 移 的 影响 


选取 拱 跨 比分 别 为 1/15、1/14、1/13、1/12、 
1/11,1/10; 而 预 张力 与 速度 差 取 No = M, = 
10 kN/m, 6(v) =40.7 m/s。 计 算出 不 同 拱 跨 比 的 
壕 形 膜 在 振动 时 产生 的 最 大 位 移 ,结果 见 表 4 和 
图 7。 


表 4 不 同 拱 跨 比 下 的 最 大 位 移 比较 
Tab.4 Comparison of the value of maximum displacement 


under different arch-span ratios 


拱 跨 比 


1/15 1/14 1/13 1⁄2 1/11 1/10 


理论 结果 /mm 9.30 9.21 9.09 8.94 8.75 8.52 


数值 结果 /mm 9.00 8.97 8.90 8.79 8.63 8.39 


相对 差 /% 3.23 2.60 2.09 1.68 1.37 1.52 


从 表 4 和 图 7 可 以 看 出 , 随 着 拱 跨 比 的 增 大 ,位 
移 逐 渐 减 小 , 拱 跨 比 从 1/15 增加 到 1/10, 膜 面 位 移 
理论 解 减少 了 8.4% , 而 数值 解 减少 了 6.8% 。 这 表 
明 拱 跨 比 与 预 张力 一 样 对 膜 面 抵抗 变形 的 能 力 有 一 
定 的 影响 ,由 于 膜 面 拱 跨 比 的 增加 ,使 膜 面 形成 了 更 
大 的 拱 状 传 力 路 径 ,为 膜 面 抵消 了 更 多 部 分 冲击 力 ， 
因此 拱 跨 比 能 提高 膜 面 承载 力 \ 增 强 膜 面 抵抗 变形 
的 性 能 。 同 时 ,理论 结果 与 数值 结果 之 间 的 相对 差 
随 拱 跨 比 的 增 大 而 减 小 ,整体 趋势 相 一 致 ,最 小 相对 


ZN 1.5296 ,最 大 为 3.23% ,符合 误差 范围 内 ,表明 


了 理论 解 的 可 行 性 。 
9.4 
一 = 一 理论 解 
一 一 数值 解 
9.0 
g 
£ 
DA 
B 8.8 


9o. 
a 


a 
> 


V5  Vi4 e Vi VH —wvio 
uis 
图 7 不 同 拱 跨 比 下 的 最 大 位 移 理论 解 与 数值 解 对 比 


Fig, 7 Comparison of theoretical and numerical solutions of 
8 


Smaximum displacement under different arch-span ratios 
513) 速度 差 对 最 大 位 移 的 影响 

忆 预 张力 以 及 摧 跨 比 取 :N= No, =10 KN/m f, = 
f S0. 1, ide 1 得 到 的 不 同 速度 差 什 与 对 应 的 模 


dr acre EHE HUC Fa lc e 2e (ERE M OE 


最 以 位 移 , 结 果 见 表 5 和 图 8。 
CN ms 不 同 速度 差 下 的 最 大 位 移 比较 


ex Tab.5 Comparison of the maximum displacement 


LE] 


under different velocity differences 


XERE2E/ (m * s7!) 


38 39 40.7 42 43 


7.96 8.16 8.52 8.79 9.00 


Seil 7mm 7.61 7.98 8.39 8.39 8.71 


相对 差 /% 4.40 2.21 1.53 0.91 0.44 


一 一 理论 解 


一 一 数值 解 


位 移 /mm 


i 39 40 41 42 E 
冲击 速度 差 /(m * s") 
图 8 不 同 速度 差 的 最 大 位 移 理论 解 与 数值 解 对 比 


Fig.8 | Comparison between theoretical solution and numerical 


solution of maximum displacement value 


for different velocity differences 
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从 表 5 和 图 8 可 以 看 出 , 随 着 速度 差 的 增加 , 产 
生 的 最 大 位 移 也 呈 增 长 的 趋势 ,这 反映 了 冲击 所 带 
来 的 能 量 越 大 , 薄膜 耗 散 该 能 量 所 需要 的 变形 也 越 
大 ;表明 了 冲击 载荷 对 薄膜 结构 影响 是 巨大 的 ,能 使 
膜 面 产生 较 大 的 变形 。 而 理论 结果 和 数值 结果 变化 
规律 一 致 ,相对 差 随 着 速度 差 的 增加 而 降低 ,最 大 差 
值 仅 4.4% 。 


6 结 it 


本 研究 利用 小 球速 度 差 模拟 得 到 冲击 载荷 , 通 
过 随机 摄 动 法 求解 得 到 了 逻 形 膜 结 构 振 动 的 位 移 函 
数 , 并 结合 算 例 分 析 与 数值 模拟 的 对 比分 析 , 得 到 以 
下 结论 。 

1) 本 人 研究 基于 随机 摄 动 法 求解 鞍 形 腊 结 构 振动 
非 线 性 微分 方程 ,可 以 得 到 理论 解 解析 ,该 理论 解 得 
能 够 预测 鞍 形 腊 在 冲击 载荷 作用 下 的 瞬 态 特性 ,得 
到 与 数值 模拟 较 吻 合 的 结果 ,说 明 本 研究 的 理论 模 
型 可 以 作为 求解 曲面 膜 结构 振动 问题 的 一 种 有 效 
方法 。 

2) 冲击 产生 的 位 移 随 预 张 力 和 拱 跨 比 的 增加 呈 
非 线性 下 降 ,但 随 着 小 球速 度 差 的 增加 而 增加 ,而 预 
张力 对 随机 振动 结果 的 影响 相对 较 大 ,这 揭示 了 鞍 
形 膜 结 构 振 动 的 瞬 态 响应 的 基本 规律 。 

3) 冲击 载荷 造成 的 薄膜 变形 是 很 明显 的 , 奏 载 
作用 带 来 的 动能 需要 薄膜 产生 较 大 的 变形 来 耗 散 ， 
严重 的 冲击 会 造成 薄膜 撕 裂 破坏 。 

膜 结 构 在 自然 环境 中 所 承受 的 风 、 暴 十 和 冰 看 
作用 具有 冲击 和 荷载 性 质 。 为 保证 膜 结 构 的 安全 性 、 
耐久 性 、 适 用 性 ,在 膜 结构 的 设计 时 ,应 设计 合适 的 
预 张力 与 拱 跨 比 ,提供 足够 抗 冲 击 变 形 能 力 。 同 时 
还 应 结合 随机 概率 统计 与 可 靠 度 理论 ,进行 膜 结 构 
在 冲击 荷载 作用 下 的 随机 动力 可 靠 度 的 进一步 研 
究 ,使 膜 结 构 的 设计 理论 更 加 合理 ,以 减少 或 防止 膜 
结构 工程 事故 的 发 生 。 
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